Spin dynamics of cold fermions with synthetic spin-orbit coupling 
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We consider spin relaxation dynamics in cold Fermi gases with a pure-gauge spin-orbit coupling 
corresponding to recent experiments. We show that such experiments can give a direct access to the 
collisional spin drag rate, and establish conditions for the observation of spin drag effects. In the 
recent experiments the dynamics is found to be mainly ballistic leading to new regimes of reversible 
spin relaxation-like processes. 
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The development of spintronics, the branch of physics 
studying spin-determined dynamical and transport phe- 
nomena was mainly related to solid-state structures. In 
these systems, spin-orbit coupling (SOC), one of the 
key elements of spintronics, is well-understood [H-Q and 
many interesting spin-related effects have been studied 
experimentally and theoretically. Very recently a de- 
tailed study of spin dynamics in ultracold atomic gases 
has become experimentally feasible In particular, 

two systems where SOC is produced by a special design of 
optical fields, attracted a great deal of attention. One of 
them is the spin-orbit coupled Bose-Einstein condensates 
[El, 0] , with the pseudospin 1 /2 degree of freedom. The 
other class is represented by the cold fermion isotopes 
40 K in Ref. Q and much lighter 6 Li studied in Ref. |{|. 
In both cases, in addition to the SOC, an effective Zee- 
man magnetic field can be produced optically. 

Typically in solids, e.g. in doped semiconductors, a dis- 
order, randomizing motion of electrons, plays the dom- 
inant role in the spin dynamics. The electron-electron 
collisions become crucial only at high temperatures, or 
in intrinsic semiconductors with optically pumped elec- 
trons and holes ■ From this point of view, cold atomic 
gases offer a unique possibility of seeing basic effects of 
interactions in the pure form since the disorder is absent 
there. The interatomic collisions lead to the spin drag 
determining the spin diffusion and, as we will see below, 
can be important for the spin dynamics in cold Fermi 
gases with SOC. 

It is well-appreciated that in the presence of strong 
SOC the effects of interatomic interactions in the spin 
dynamics are difficult to analyze as this requires tracing 
essentially coupled orbital and spin subsystems. Fortu- 
nately, these dynamics become uncoupled not only for 
vanishing SOC, but also when it corresponds to an ef- 
fective non-Abelian vector potential 1T23J of a pure 
gauge form, which happens in a broad class of systems. 
Remarkably, the three-dimensional (3D) fermionic gases 



with SOC realized in recent experiments @,Sli3| belong 
to this interesting class. For a pure gauge SOC the be- 
havior of the physical system maps to that of a system 
without SOC, which allows to consider effects of SOC 
of an arbitrary strength. In this case all qualitative fea- 
tures of the spin dynamics are the same as for a generic 
SO field, but the analysis is much easier. Here we study 
spin dynamics for systems with a pure gauge SO cou- 
pling, where the entire pattern even if it is complicated 
by the interatomic interactions, can be explored using a 
single formula. As a byproduct of our analysis we show 
that studying the spin relaxation in SO coupled systems 
may offer an alternative route to experimentally access 
the spin drag resistivity in SOC-free systems. 

We represent the Hamiltonian of a 3D Fermi gas of 
atoms with the mass m and a linear in momentum SOC 
as follows 
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where ip is a spinor fermionic field operator, a a with 
a = {x, y, z} are the Pauli matrices, the space com- 
ponents Aj of an effective SU(2) potential parametrize 
SOC, and B b is the Zeeman field (the time component 
of the SU(2) potential) (H- The term H int [ip\ ip] in 
Eq. ([T]) describes the interparticle interactions and the 
external trapping potential. In the following we assume a 
short-range interaction characterized by a s-wave scatter- 
ing length a s . For simplicity we consider a homogeneous 
system, assuming the validity of the local approximation. 

The spin relaxation is probed by applying the Zeeman 
field of the form B(t) = B0(— t). That is, we first bring 
the system to an equilibrium polarized state in the pres- 
ence of the static Zeeman field which is then switched off 
at t = 0. After the polarizing field is released, the spin 
polarization relaxes to zero because of the combined ac- 
tion of the SOC-induced spin precession and the particle- 
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particle collisions. 

In conventional semiconductors the SOC is usually 
weak, that is | A°j | -C pf , where pf is the Fermi momen- 
tum. As a result, only the vicinity of the Fermi surface 
is important for the spin dynamics. Moreover, the weak 
coupling practically does not influence the initial spin- 
polarized state. With optical fields one can achieve the 
ultrastrong coupling where the effective SU(2) field \ A®\ 
is close to or even larger than the Fermi- or the mean ther- 
mal momentum. In this ultrastrong coupling regime the 
collisional relaxation processes can become irrelevant and 
the dynamics is mainly due to the nonuniform spin pre- 
cession, that is its broadening in the momentum space. 
In addition, the strong SOC significantly modifies the 
initial state making it strongly non-semiclassical. 

While the description of the spin dynamics for a strong 
generic SOC is very complicated, significant simplifica- 
tions occur in the special case of a pure gauge SO field 
of the following "direct product" form 



f 



M = -q^ 
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Here n is a unit vector pointing along the SO precession 
axis (momentum independent in this case), and the vec- 
tor q determines the strength of SOC for every direction 
of the particles' momentum coupled to the spin degrees 
of freedom. Exactly this type of the SU(2) field has been 
realized in the recent experiments with cold Bose and 
Fermi gases [ISH- 

Because different components of the SU(2) potential in 
Eq. ([5]) commute, it can be gauged away by a local SU(2) 
rotation of the fermion fields: 



■0(r) = ^( r )'0( r ) 5 U(r) = exp 
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After the rotation, the Hamiltonian of Eq. ([I) reads 
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where B(r) = tr{oi/ _1 (r)(erB)W(r)}/2 is the trans- 
formed Zceman field. For definiteness we consider the 
configuration of fields generated in the experiments of 
Ref. pi . This corresponds to n = (0,0, 1), q = (q, 0,0), 
and B = (B, 0, 0), which translates to the following heli- 
coidal structure 261-291 of the transformed Zeeman field 



B(r) = B [xcos(qx) +ysin(gx)] 



(5) 



Since the interparticle interaction and the trap potential 
are spin independent the form of H ln t is untouched by the 
SU(2) rotation of Eq. Q. Therefore the physical system 
with the pure gauge SOC subjected to a uniform polar- 
izing field B = x_B is mapped to a dual system without 
SOC in the presence of an inhomogeneous Zeeman field 
- the magnetic helix of Eq. ([5J . 



The mapping of the spin relaxation for a pure gauge 
SOC to the diffusive spin dynamics in the SO-free sys- 
tem has been already noted in the context of conven- 
tional semiconductors 22, [23| where the spin relax- 



ation/diffusion is mostly determined by a disorder. Ap- 
parently this is an absolutely general statement, which 
in particular implies that the spin diffusion coefficient 
D s and the spin drag rate r s in cold atomic gases can be 
measured by monitoring the dynamics of a global aver- 
aged spin, provided the proper SOC in created. 

After the transformation of Eq. ([3]) the spin dynam- 
ics can be easily found using the linear response theory 
for homogeneous, SOC-free systems. Namely, we con- 
sider a response to a week Zeeman field Eq. ([5]) with 
B(t) = B9(—i). The spin density S(r, t) in the trans- 
formed system follows the helicoidal pattern of Eq. ([5]) 
with a time-dependent amplitude S(t): 



S(r, t) = S(t) [xcos(gx) +ysin(ga;)] 



(6) 



The amplitude S(t) of the spin helix can be related to the 
standard spin response function x <T<T (g,w) at the wave 
vector q (see, e.g., Ref. [30l |) 
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To recover the spin dynamics in the original system we 
apply the inverse SU(2) transformation to the spin dis- 
tribution S(r, t) in Eq. This unwinds the helix and 
leads to the evolutions of the spatially uniform physical 
spin distribution S(t) = 5(f)x presented below. 

Equation (J7J| relates spin dynamics in the physical sys- 
tem with SOC to the smearing out of an inhomogeneous 
spin texture in a system without it. We emphasize that 
Eq. ([7]) is exact and valid for any phase without spon- 
taneous ferromagnetic ordering, either normal or super- 
fluid one. The SOC enters the spin dynamics solely via 
the wave vector dependence of the spin-spin correlation 
function x<ya(q,^) of the SOC-free system. Hence by 
monitoring the spin evolution for different values of q 
one can access the entire momentum and frequency de- 
pendence Of Xcrcr(q,w)- 

We begin with an illustration of the simplest universal 
type of spin dynamics occurring for small SO couplings, 
w a t q <C T s , where u a t is the typical atomic velocity such 
as the Fermi velocity vf or the mean thermal velocity 
vt in the degenerate and nondegenerate system, respec- 
tively. In this case the expression in the square brackets 
in Eq. ([7]) takes the diffusive form iu>/{iu> + D s q 2 ), which 
is guaranteed by the spin conservation law. Hence the 
spin relaxes exponentially S(t) = 5(0) exp (—D s q 2 t) by 
the Dyakonov-Perel' mechanism (3f|. The relaxation rate 
for a given SOC is directly determined by the spin diffu- 
sion coefficient of the dual SOC-free system. 

For an arbitrary strong coupling the dynamics is much 
richer. To explore the behavior in the full range of SOC, 



3 



densities, and temperatures, we employ a spin version 
of the "conserving" relaxation time scheme by Mermin 
[3^ | . In this scheme one approximates the collision inte- 
gral in the linearized kinetic equation for the spin part 
s p (q, t) of the Wigner function by the following form, 
Icoi = -r 8 [s p (q,t) - s p q (q,t)], where s p q (q,<) is the 
equilibrium function corresponding to the instantaneous 
value of the average spin S(q, t). For any inhomogene- 
ity wave vector q this approximation recovers the correct 
static and high-frequency response, and respects the local 
spin conservation. The dynamical spin response function 
in this scheme takes the form 



Xo(q,uJ + iT s ) 



Xo(<7,0) 



(8) 



where xoiq^) is the Lindhard function 

/p+q/2 - /p-q/2 



Yo 



(3,2) =£ 



U ~ ( £ P+q/2 - £ p-q/2) 



(9) 



with / p being the Fermi function, and e p = p 2 /2m. 

The spin drag relaxation rate T s enters the Mer- 
min's scheme as a phenomenological parameter originat- 
ing from interatomic collisions. Although collisions con- 
serve the total momentum, they cause a friction between 
different spin species, and thus lead to relaxation and 
damped response in the spin channel. Strictly speaking, 
the spin drag rate r s depends on q. It is however clear 
that for large q > p F the average spin relaxes mostly 
because of fast nonuniform precession (ballistic motion). 
The spin drag contribution becomes qualitatively impor- 
tant in the opposite limit of q <C pp. Therefore below 
we adopt the standard expression based on the Born ap- 
proximation at q = [33, 34 1 



r s = 



k 2 



^— ' 3mnT 
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dm [Imxo(fc,^)] 2 
7T sinh 2 (w/2T) 
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where T is the temperature, and n is the concentration 
of the atoms. The spin drag rate of Eq. (fl0|) has the 
following characteristic temperature dependence 



T s ~ E F (p F a s ) x 



(T / Ey) , 

(T/E F ) 1/2 



T <C E F 



T^>E F . 



(11) 



The low-T behavior is typical for the Fermi liquid, while 
in the high-T limit it describes the frequency of collisions 
between particles with a scattering cross-section ~ a 2 
moving at the mean thermal velocity vt ~ y/T/m. For 
the repulsive interaction the Stoner instability criterion 
limits the product p F a s < ir/2 j35j . thus preventing one 
from having a fast relaxation in the degenerate gas. 

Equations (|7|)- (fT0| provide us with a simple unified 
description of the spin dynamics in a wide range of tem- 
peratures, SO couplings, densities and interactions. This 



scheme captures all main physical effects which influence 
the spin relaxation. The difference in the Fermi functions 
in Eq. ^ describes the modification of the initial state 
by the SO coupling. It selects a part of the momentum 
space with the particles polarized by the initial Zeeman 
field, as it is shown in Fig. [1] for the very strong coupling. 
The difference in energies in the denominator in Eq. (|9]) 
determines the frequency of the spin precession, while the 
imaginary shift of u by iT s describes the spin drag effect. 
Now we use Eqs. ([7|)- (fT0|) to explore different regimes of 
spin dynamics and discuss those not achievable in the 
solid-state experiments. 
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FIG. 1: (Color online.) Schematic plot of precession of spins 
of atoms with different momenta s p (t), as shown by short 
arrows, in the case of very strong SOC with q > 2pF- Filled 
circles show spin-split thermally broadened Fermi spheres. 

We consider first the low-temperature degenerate gas, 
T <C E F - If SOC is weak, so that q <C pf, only particles 
near the Fermi surface participate in the dynamics. The 
characteristic time scale is set up by l/qv F , that is the 
time for ballistically traversing the period of the spin he- 
lix. The relation of this time to the spin drag rate T s de- 
termines the regime of the spin dynamics. If qv F /T s 3> 1 , 
the dynamics is ballistic and the evolution of the total 
spin is dominated by the inhomogeneous p-dependent 
spin precession of the Fcrmi-surfacc particles. Appar- 
ently this is always the case at sufficiently low tempera- 
tures, T — > 0. In the opposite limit of qv F /T s <C 1 the 
dynamics is diffusive, which corresponds to the conven- 



tional Dyakonov-Perel' spin relaxation mechanism [31 
From Eq. (|11|) we find that with the increase of the tem- 
perature the ballistic regime crosses over to the diffusive 
one at Tdb ~ {v F / a s )\J q/PF- Condition T^b <C E F im- 
plies q <C p F (p F a s ) 2 for the SOC strength. Taking into 
account that, e.g., in Ref. 0, PFa s ~ 10 -2 , we conclude 
that to reach the diffusive regime one has to increase the 
interaction or to make the SOC considerably weaker. 

Apart from determining the spin drag rate, the tem- 
perature effects modify the dynamics via broadening the 
Fermi distribution in Eq. (|9]). Indeed, the broadening of 
the velocity distribution by Sv F ~ (T/E F )v F causes the 
broadening in the precession rate by the typical value of 
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FIG. 2: (Color online.) Time dependence of the total spin at 
strong spin-orbit coupling, (a) low temperature, degenerate 
gas, (b) high temperature, classical gas. The values of q and 
T are shown in the plots. 



SD, = (T/Ep)qvp. In the ballistic regime the resulting 
time dependence of the total spin is 



S(t) = S(0) 



sinh(7rf(5ft) Ep ' 



(12) 



where the rapid increase in si nh( 7rf(5fj) at t > 1/SCl leads 
to the exponential damping |36| in the spin oscillations 
initially following the sin(vpgt) / vpqt time dependence. 
At Sfl/r s < 1, that is at T > T% b /E F the asymptotic 
behavior is determined by T s rather than by <5f2. 

For strong SOC the dynamics is mostly ballistic, but 
still quite nontrivial and diverse. As it can bee seen in 
FigUJa), at q = pf the spin relaxes to zero monotonically, 
while pronounced large amplitude oscillations appear for 
ultrastrong SOC with q^ pp. The reason is that in the 
latter case we have two well separated spin-split Fermi- 
spheres. The position of the sphere center determines the 
mean precession rate, and the broadening is determined 
by the Fermi momentum pp. An interesting intermedi- 
ate regime occurs at q = 2pp when the spin-split Fermi- 
surfaces touch each other and the time dependence of the 
spin takes the form 



S(t) = S(0) 



sin 2 (4E F i) 
(4E F t) 2 ' 



(13) 



on the temperature [see Eq. (fTTj) ]. Therefore the dynam- 
ical regime is temperature-independent - the dynamics is 
ballistic if (q/pp)/ (PF a s) 2 ~S> 1 and diffusive otherwise. 
The condition of diffusive relaxation is equivalent to the 
condition Tdb -C Ep. As a result, if q > pp (ppa s ) 2 , the 
spin dynamics is ballistic at any temperature. In this case 
the width of the thermal Maxwell distribution is much 
larger than the scale of SOC and the inhomogeneous pre- 
cession is equivalent to the Fourier transform of the Gaus- 
sian momentum distribution with the coordinate qt/m. 
For a weak SOC in the ballistic regime the spin relax- 



(v T qt) /6 



with 



Now wc consider the high-temperature case, T 3> Ep. 
In this limit the mean free path ~ vt /r s does not depend 



ation is Gaussian S(t) = 5(0) cxp 

the time scale of 1/qvr- At larger scattering length a s , 
that is at larger T s , we enter the diffusive regime where 
the Gaussian damping slows down to the exponential 
Dyakonov-Pcrel' relaxation with the timescalc T s /q 2 v^. 
In the nondegenerate gas the function S(t) is practically 
always monotonic. To detect signatures of the oscillatory 
spin precession one needs to go to the ultrastrong SOC 
and to satisfy experimentally the condition <}> pt- 

It is worth noting that in general for the weak SOC 
the spin dynamics in the ballistic limit is determined by 
the parameter v ai q since this is the only characteristics 
of the low-energy excitation spectrum entering Eq. (|9j). 
For the strong coupling with clearly separated spin-split 
momentum distributions we have two main parameters 
describing spin dynamics: spin precession rate of the or- 
der of q 2 /2m and spin relaxation rate of the order of qv a t, 
as can be seen from comparison of Figl2ta) and (b) for 
the low- and high-temperature gases, respectively. This 
links the spin relaxation dynamics to the spectrum of 
electron-hole excitations in the Fermi gas. 

We emphasize that the above limiting regimes of 
the spin dynamics are described within a single uni- 
fied scheme, which is completely determined by four 
experimentally controllable parameters, the density n, 
the temperature T, the scattering length a s , and the 
SOC strength q. Apparently all crossovers and transi- 
tion regimes are also covered by this approach. 

To conclude, we have shown that the spin dynamics 
in cold Fermi gases can be mapped on the spectrum of 
particle-hole excitations. The spin drag rate can be ob- 
tained from the relaxation of the uniform spin density 
without generating an inhomogeneous spin distribution 
in the real-space, provided the conditions of diffusive spin 
dynamics in terms of the SOC strength, concentration, 
and the scattering length are satisfied. 
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